Condition of the occurrence of phase shp centers in superconducting nanowires under 

apphed current or voltage 
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Experimental results on the phase slip process in superconducting lead nanowires are presented 
under two different experimental conditions: constant applied current or constant voltage. Based 
on these experiments we established a simple model which gives us the condition of the appearance 
of phase slip centers in a quasi-one-dimensional wire. It turns out that the competition between two 
relaxations times (relaxation time of the absolute value of the order parameter t\^\ and relaxation 
time of the phase of the order parameter in the phase slip center r^) governs the phase slip process. 
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I. INTRODUCTION 

Since the discovery of superconductivity it was ex- 
pected that if the superconductor was subjected to a 
constant electric field superconductivity will inevitably 
be destroyed. The reason is that the superconducting 
electrons will be accelerated by the electric field and will 
reach a velocity above its critical velocity. However if the 
sample is short enough or if the electric field exists only in 
a small part of the sample such that the path over which 
the Cooper pairs are accelerated are sufficiently short an 
electric field can exist in the sample in the presence of 
superconductivity. Another example is the presence of 
an electric field in the superconducting sample which is 
attached to a normal metal. In this geometry the in- 
jected current from the normal metal will be converted 
into a superconducting current on a distance of about 
the charge imbalance distance (Aq) and on that scale 
an electric field will exist in the samplei. In this case 
the electric field is compensated by the gradient of the 
chemical potential /ig of superconducting electrons and it 
does not lead to an acceleration of the superconducting 
condensate (see for example the book of Schmidt^). 

But there is another mechanism which allows super- 
conductivity to survive in the presence of an electric field 
deep in the superconducting sample of arbitrary length. 
This is the phase slip mechanism. Initially this phenom- 
ena was used in order to estimate the relaxation time of 
superconducting current in a superconducting wire?. If 
the order parameter vanishes in one point of the wire the 
phase of the superconducting order parameter exhibits a 
jump of 2n at that point- and as a result the momentum 
p ~ decreases by 2-k/L. So even if the electric field 
accelerates the superconducting electrons it does not lead 



to a destruction of supercondictivity because the momen- 
tum is able to relax through the phase slip mechanism. 

This simple idea was understood already a long time 
ago. Therefore, it is very surprising that there are prac- 
tically no experimental nor theoretical works studying in 
detail what will happen if a voltage (i.e. electric field) is 
applied to the superconductor. In previous works mainly 
the situation with applied current (i.e. the I=const 
regime) was studied. In the latter case the phase slip 
process was studied theoretically in detail (see for ex- 
ample the review of Ivlev and Kopnin'* and the books of 
Tinkham.^ and Tidecks5) . On the basis of a numerical so- 
lution of the extended time-dependent Ginzburg-Landau 
equations it was found that the phase slip (PS) phenom- 
ena exists in some region of currents. The lowest critical 
current at which a PS solution first appears in the sys- 
tem may be smaller than the depairing Ginzurg-Landau 
current density and as a result it leads to hysteresis in 
the I-V characteristics in the I=const regime^. 

The I=const regime was also studied experimentally in 
a number of works'*'^ starting from the paper of Meyer 
and MinnigerodaSii. The most characteristic effect of 
the phase slip mechanism is the appearance of a stair 
like structure in the current-voltage characteristics. In 
Ref.^ a simple phenomenological model was built in order 
to quantitatively describe this feature. It was proposed 
that every step in the I-V characteristic is connected with 
the appearance of a new phase sip center (PSC) that in- 
creases the resistivity of the sample by a finite value. In 
this model, it was conjectured that in a region of the 
sample with size of about the coherence length ^ fast os- 
cillations of the order parameter occurs in time which 
produces normal quasi-particles. Because of the finite 
time needed to convert normal electrons into supercon- 
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ducting electronsiiS there is a region of size of about Aq 
near the phase slip center where the electric field and 
the normal current are different from zero. It means that 
the chemical potential of superconducting fi^ and normal 
electrons /i„ are different near PSC and their difference is 
proportional to the charge imbalance Q in a given point 
of the superconductor^-*. Because usually Aq ^ it 
is possible to neglect the ^-sized region with oscillating 
\ip\ and consider that phenomena as a time- independent 
process. Such a model gives a finite voltage, which is 
connected with every phase slip center, which is equal 
to Vpsc — '^A.QPn{I — pic)/S with pn the normal resis- 
tivity, S the cross-section of the current carrying region, 
Ic is the critical current and /3 < I is a phenomenologi- 
cal parameter. In the experiment of Dolan and Jackel^ 
the distribution of fj,s and fin near a PSC were measured 
which fully supported the idea of Ref.? that a difference 
exist between /is and /i„ near the phase slip center. 

Despite numerous theoretical and experimental works 
the physical conditions under which phase slip phenom- 
ena can exist is still not clear. In a review'* on this subject 
it was claimed that phase slip phenomena are connected 
with the presence of a limiting cycle in the system which 
leads to such type of oscillations. But it was not ex- 
plained how and why it leads to the phase slip phenom- 
ena. 

Based on our previous work-'^Si on the time-dependent 
Ginzburg-Landau equations (the investigated systems 
were superconducting rings in the presence of an external 
magnetic field) we know that systems which are governed 
by such equations exhibit two relaxation times. One is 
the relaxation time of the phase of the order parameter 
and the other is the relaxation time of the absolute 
value of the order parameter T|^|. In Rcf.*" it was estab- 
lished that phase slip processes can occur in such systems 
when roughly speaking < T|^|. In the present paper 
we discuss this question in the context of superconduct- 
ing wires in the presence of an applied current or voltage. 

To our knowledge there exists only a single theoretical 
work in which a superconducting wire in the presence of 
an applied voltage was studiedii. The authors used the 
simple time-dependent Ginzburg-Landau equations and 
found that the behavior of the system is very compli- 
cated and strongly depends on the length of the wire and 
the applied voltage. Neither a detailed analysis nor any 
physical interpretation of their results was presented. In 
a recent letter^^ we presented our preliminary theoretical 
and experimental results on the dynamics of the super- 
conducting condensate in wires under an applied voltage. 
It turned out that in this case the I-V characteristics ex- 
hibit a S-shape which was explained by the appearance 
of phase slip centers in the wire and their rearrangement 
in time. In the present paper we will present more de- 
tails and extend our previous work to the situation in 
which defects are present, and we investigate the effects 
of boundary conditions and an applied magnetic field. 

The paper is organized as following. In Sec. II we 
present our theoretical results and give the conditions for 



the existence of phase slip centers when current (Sec. Ila) 
or voltage (Sec. lib) is appHed to the superconducting 
wire. In Sec. Ill we show our experimental results and 
in Sec. IV we compare theory and experiment. 

II. THEORY 

We study the current-voltage characteristics of quasi- 
one-dimensional superconductors using the generalized 
time-dependent Ginzburg-Landau (TDGL) equation. 
The latter was first written down in the work of Ref^ 

= (V-iA)V + (l-|^n^. (1) 

In comparison with the simple time-dependent Ginzburg- 
Landau equation where 7 = it allows us to describe a 
wider current region (with a proper choice of the pa- 
rameters u and 7) where a superconducting resistive 
state exists and gives us a wider temperature region in 
which Eq. (I) is applicablai2*i^*i^. The inelastic col- 
lision time te for electron-phonon scattering is taking 
into account in the above equation through the tem- 
perature dependent parameter 7 = 2TE^o(T)/h (Aq = 
iksTcU^^^il - T/TcY/^/t: is the equihbrium value of the 
order parameter). 

Eq. (1) should be supplemented with the equation for 
the electrostatic potential 

Aip = div(Im(i/)*(V-iA)V')), (2) 

which is nothing else than the condition for the conser- 
vation of the total current in the wire, i.e. divj — 0. 
In Eqs. (1,2) all the physical quantities (order param- 
eter ip = iV'le*'^, vector potential A and electrostatical 
potential ip) are measured in dimensionless units: the 
vector potential A and momentum of superconducting 
condensate p ~ Wcj) — A is scaled by the unit $o/(27r^) 
(where $0 is the quantum of magnetic flux), the or- 
der parameter is in units of Aq and the coordinates 
are in units of the coherence length £,{T). In these 
units the magnetic field is scaled by Hc2 and the cur- 
rent density by jo = c^q/Stt^A^^. Time is scaled in 
units of the Ginzburg-Landau relaxation time tgl = 
47r(T„A^/c^ — 2Th/nAl, the electrostatic potential ((/?), 
is in units of (po — c<I>o/87r^^A(j„ = h/2eTGL {o'n is 
the normal-state conductivity). In our calculations we 
mainly made use of the bridge geometry boundary con- 
ditions: \^{-L/2)\ = |^(£/2)| = 1, ip{-L/2) = 0, 
ip{L/2) = V and ^lj{L/2,t + dt) = V'(i/2, t)e-*'^(-^/2)dt^ 
We chose these boundary conditions because at low tem- 
peratures the normal current is converted to a supercon- 
ducting one due to the Andreev reflection on a distance of 
about ^0 — O-lShvp/kBTc near the S — N boundarjii2i. 
It means that practically there is no injection of quasi- 
particles from the normal material to the superconductor 
and hence we can neglect the effect of charge imbalance 
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near the S-N boundary. The bridge geometry bound- 
ary conditions models this situation. The parameter u is 
about 5.79 according to Refpi^. We also put A — Qin Eq. 
(1,2) because we considered the one-dimensional model 
in which the effect of the self-induced magnetic field is 
negligible and we assume that no external magnetic field 
is applied. 



A. Constant current regime 

Lets consider first the more simple case when a con- 
stant external current is applied to the sample. In such 
a case it was theoretically foun d^^i^^ that the system ex- 
hibits hysteretic behavior. If one starts from the super- 
conducting state and increases the current the supercon- 
ducting state switches to the resistive superconducting 
or normal state at the upper critical current density jc2 
which, in a defectless sample, is equal to the Ginzburg- 
Landau depairing current density jc l = v'4/27jo ■ Start- 
ing from the resistive state and decreasing the current it 
is possible to keep the sample in the resistive state even 
for currents up to jd < jc2 (which we call the low critical 
current). For jd < j < jc2 such a state is realized as a 
periodic oscillation of the order parameter in time at one 
point of the superconductoriiii^. When the order param- 
eter reaches zero in this point a phase slip of 27r occurs. 
This is the reason why such a state is now called a phase 
slip state and this point a phase slip center (PSC). Using 
results obtained in our earlier work^" we claim that the 
value of jci depends on the ratio between the two char- 
acteristic times in the sample: the phase relaxation time 
of the order parameter and the relaxation time of the 
absolute value of the order parameter in the region 
(with size of about ^) where the oscillations of the order 
parameter occurs. 

Having written Eq. (1) for the dynamics of the phase 
and the absolute value of the order parameter in a quasi- 
one dimensional wire of length L (-L/2<s<L/2) 

uVTT^^ = ^ + mi-m'-p% (3a) 



dt 



LP - 



ds 



(3b) 



it is easy to estimate both these times. Indeed, from Eq. 
(3a) it directly follows that 



-I- 7^|'0p — wy (for7 ^ 1). 



(4) 



To determine we need to know how fast the phase 
changes over the region (with size of about ^) at which 
the order parameter oscillates. Our numerical calcu- 
lation shows that the time-average of the derivative is 
d{(l))/dt = Ce{s) (with C(s)=const and e{s) the theta 
function). Because the time-averaged electrostatic po- 
tential ip changes over a distance Aq ^ 1 for 7 3> 1 
(where Ag = yjl -f 72|i/)|2/u|i/)|2 ~ -y/y jg tJ^e decay 



length of the normal current density and the charge im- 
balance Q) and that the order parameter is about unity 
already at s = we can estimate this constant as 
C ~ Aqdjn/ds ~ j I-^q- Consequently we find that 



r0 



Agj 



(5) 



The product Agj is roughly the voltage drop Vpsc 'pro- 
duced' by the phase slip centes^. Thus the time change 
of the phase of the order parameter in the phase slip 
center is inversely proportional to the voltage drop over 
the whole structure (and C = Vpsc/2). This is a 
consequence of the fact that the time-averaged electro- 
chemical potential of the superconducting electrons /zf = 
d{(j)) /dt does not change in space and its jump near the 
phase slip center is equal to the voltage drop of this struc- 
ture in the sample. In terms of the language used by 
Schmid and G. Schorki^, t\^\ is the 'longitudinal' time 
and decreases with increasing 'transverse' time be- 
cause ~ 1/Aq ~ \l ^Jtq. In RefiiS it was found 
that the phase slip events occur periodically in time when 
< T|^|. This gives us an estimation for the lower crit- 
ical current 
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Jcl 



"I^|Aq 



(6) 



By varying the parameter 7 we can change both t\^\ 
and Aq and hence we can vary jci- In Fig. 1 we show 
the voltage in the sample, the normal current density and 
the absolute value of the order parameter in the phase 
slip center averaged over the time as a function of the 
external current for two different values of 7. In our sim- 
ulations we started from the superconducting state with 
j < 3gl- At j > jgl the system jumps instantaneous to 
the resistive state with finite voltagai^. Than we decrease 
the external current and at j < jci the system transits 
to the purely superconducting state with y = 0. We 
should emphasize that at j — > jci the voltage jumps by a 
finite value (this result is qualitatively different to the re- 
sults of Ref>i2*i^ where the authors found that y ^ at 
j jcl-)- It means that there is a non-infinite maximal 
oscillation period for the order parameter in the phase 
slip center. We believe that the finite voltage jump AV 
or finite period of oscillations is directly connected with 
the threshold condition ~ T|^| for activation of regular 
phase slip processes in the constant current regime and 
it means that AV ~ 1/t|^|. 

Before going further we should stress here that the 
above condition for the existence of a phase slip process 
T0 < r|^| is a rather rough estimate. Indeed, Eqs. 3(a,b) 
are a coupled system of equations, and besides r|^,|, it 
exphcitly (see Eq. (4)) depends on the value of How- 
ever, the above condition allows us to explain the general 
qualitative properties of the phase slip processes (includ- 
ing the existence of AV, jd and its dependence on 7 
and u) and predict new features which will be discussed 
below. 
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FIG. 1: The dependence of the time averaged voltage (a) 
on the external current for a wire containing only one phase 
slip center. In figures (b,c) the dependencies of the normal 
current density (b) and the order parameter (c) in the phase 
slip center are shown. Length of the wire is 40^. Solid curves 
correspond to 7 = 10(Aq ~ 2.3^), and the dotted curves are 
for 7 = 40 (Aq ~ 4.10. 



Far enough from jd the dependence of V{j) on the 
current is close to linear. When the current increases, 

decreases (see Eq. (5)) and hence the time Tpsc be- 
tween two phase slips also decreases and the order param- 
eter has less time for recovering at the phase slip center. 
That is the reason why the time averaged voltage in- 
creases (as {V) = 2n/Tpsc), the averaged order parame- 
ter (I '0(0) I) decreases and the fraction of the normal cur- 
rent {jn{0))/j increases with increasing external current. 
It is interesting to note that from the phenomcnological 
Skocpol-Beasley-Tinkham (SBT)? model it follows that 
(j„(0))/j = 1 — Pjci/j which quahtatively resembles the 
dependence shown in Fig. 1(b). 

Let us now discuss the effect of defects on the I-V char- 
acteristic. This question was considered previously in 
Ref. for two different models of defects: local variation 
of the critical temperature and a local variation of the 
mean free path. We will repeat these calculations par- 



tially and interpret it in terms of a competition between 
and T|^|. In addition to the first type of defect we will 
also study the effect of the variation of the cross-section 
of the wire. 

The first type of defect is the inclusion of a region 
in the superconductor which suppresses Tc and the or- 
der parameter becomes lower than the equilibrium value 
Aq even in the absence of any external current. This is 
modelledSS by introducing the term p{s)ijj in the right 
hand side (RHS) of Eq. (1). In the present calculation 
we choose p{s) — —po0{O.5 — \s\). The larger po the 
more the order parameter is suppressed in the center of 
the wire. Firstly, such a defect leads to a decrease of 
the upper critical current jc2- Secondly, it decreases the 
lower critical current jd- Indeed, when we introduce a 
defect, the RHS of Eq. (3a) decreases at the spatial posi- 
tion where the order parameter oscillations and hence \ip\ 
needs more time to change in that spatial position. Thus 
this type of defect leads to an increase of the relaxation 
time of the order parameter in the region of the defect 
(as an indirect prove of this we found a decreasing AV^ 
with increasing 'strength' of the defect). If the size of 
the defect is smaller than Aq we can neglect its effect on 
the relaxation time of the phase of the order parameter. 
Finally, we can conclude that the 'stronger' the defect 
the smaller the value of the lower critical current jci ■ 

In Fig. 2(a) we present our numerical results for two 
different values of pg. With increasing 'strength' of the 
defect the upper and lower critical currents decrease and 
start to merge. As a consequence the hysteresis in the 
I-V characteristics will disappear when the defect is suf- 
ficiently strong. 

The second type of defect is one for which we have a 
local decrease of the cross-section of the wire. In this 
case we expect that r|^| will not be influenced. The sit- 
uation with T0 is more complicated because a large part 
of the superconductor with size of about Aq participates 
in the formation of this time. Two limiting cases can 
be distinguished. If the defect is much smaller than Aq 
we can neglect its effect on and hence we will have 
the same lower critical current as for the case of an ideal 
wire. In the opposite case of a defect with size much 
larger than Aq we only have to take into account the 
increased current density in the part of the wire where 
we have a smaller cross-section and apply Eq. (5). In 
this case the current jd is decreased by a factor Dav/Di, 
where Dav is the average cross-section and Di < Dav is 
the local reduced cross-section. 

Another important question is the value of the upper 
critical current jc2- If the size of the defect is larger 
than ^ then the proximity effect from adjacent parts near 
the defect will be small and jc2 will decrease by a factor 
Dav/ Di (in the opposite case the current jc2 is almost 
defect independent). Therefore, for a defect with length 
^ < ^ ^ Aq the I-V characteristic may be reversible with 
a proper choice of the parameters (like in the case of a 
local suppression of Tc). In Fig. 2(b) we show the results 
of our numerical calculations and we obtain a qualitative 
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to take a smaller time step (see Eq. (1)) which increases 
the computation time. 
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FIG. 2: The dependence of the time averaged voltage on the 
external current for a wire containing a single defect. Figure 
(a) corresponds to a wire with a local suppression of the Tc 
(dotted curve for po = —2, dashed curve for p = and solid 
curve for a wire without defect). Figure (b) corresponds to 
a local variation of the cross-section of the wire which we 
modelled as D{s) = 1 - fle"""/'^. Solid curve in this figure 
is for a wire without a defect, dashed curve for a wire with 
defect parameters /3 = 0.5 and Id = 0.5, dotted curve with 
f3 = 0.2 and Id — 2, dash-dotted curve with (3 = 0.5 and 
Id = 2. The length of the wire is L = 40^ and the parameter 
7=10 (Aq~2.30. 



agreement with the above physical arguments. When the 
length of the defect is smaller than Aq then the lowest 
critical current density practically does not change. In 
the opposite case jd decreases by a factor of Dav/Di. 
The upper critical current density changes considerably 
only if the length of the defect exceeds ^. Unfortunately, 
we are not able to consider the case for which the length 
of the defect is simultaneously much smaller than Aq and 
much larger than ^ because of computational restrictions. 
For example, when we increase 7 by a factor of two the 
time of calculations increases also by a factor of two but 
the ratio Ag/^ increases roughly only by a factor of \/2 
for 7 3> 1 . The reason is that when we increase 7 we need 
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FIG. 3: Dependence of the low critical current jd on the 
length of the superconducting wire for different values of 7. 
The current jci is normalized to its value (.7*1) at lengths 
L » Aq. With increasing 7 the decay length of the normal 
current increases and hence jd starts to depend on L for 
longer wires. In the inset, the dependence of jci on the length 
of the superconducting wire for S-N boundary conditions is 
shown. 

When the wire has a length less than Aq then this 
will affect the distribution of the normal current density 
in the wire and hence the time and the lower critical 
current jd ■ Indeed from our bridge boundary conditions 
it follows that djn/ds{±L/2) = (see Eq. 3(b)). Then 
taking into account that d{(j)) / dt = V9{s) /2we can easily 
obtain, in the limit Aq ^, that 



On(0)) 



1 



2Aq tanh(L/2AQ) 



(7) 



where the current (jn(0)) is the average normal current 
in the phase slip center. 

The voltage jump ~ 2'k/t\^\ at jd should not 
change with varying L (at least for L ^ when prox- 
imity effect does not have an effect on t\^\)- From Eq. 
(7) it follows then that the normal current density in- 
creases with increasing length. But (in(0)) should al- 
ways be smaller than the full current j. It implies that 
the critical current density jd will increase with increas- 
ing length of the wire as jd ^ l/tanh(L/2AQ) in order 
to keep the ratio (j„(0))/j constant. Fig. 3 illustrates 
the dependence of jd on the length of the wire which 
we obtained on the basis of a numerical solution of Eqs. 
(1,2). Unfortunately, in our calculations we are not able 
to use very large values of 7 and the maximal value of 
Aq was 4.1^ for 7 = 40. But nevertheless, we found that 
jcA increases with decreasing wire length. We should add 
here that we also found that AF also decreases a little 
bit. We connect it with a small change in T|^| due to the 
small length of the wire and hence the increased effect of 
the boundaries. 
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The more pronounced effect of the finite length of 
the wire on the value of jci is found for the case that 
we used N-S boundary conditions: ip{ztL/2) = and 
d(p/ds{±L/2) = —j. These boundary conditions are ap- 
proximately valid for samples at temperatures close to 
Tc. It is easy to show that in this case 



On(0)) = 



V{j) tanh(£/2AQ) 
2Aq 1 - 1/a • cosh(i/2AQ) ' 



(8) 



where a = {jn{0))/j < 1 (see Fig. 1(b)). In this case jci 
also increases with decreasing wire length (see inset in 
Fig. 3) . But besides there is a finite length Lq for which 
(j„(0)) ^ oo. It implies that the phase slip process is 
not possible in wires with length L < Lq. In such wires 
the system goes from the superconducting state directly 
to the normal state. 

Not only the finite length of the sample is able to 
change the value of jd- If we apply a magnetic field 
parallel to the length of the wire it will suppress the or- 
der parameter in the sample. Our analysis shows that if 
the diameter of the wire d is less than 2^ and if we can 
neglect screening effects (A > ^) then the distribution 
of the order parameter will be uniform along the cross- 
section of the wire. The order parameter depends on H 
as 

IVP - 1 - {H/H,f 

with He — 2.9^Q/TT^d. This behavior is very similar 
to the behavior of a thin plate in a parallel magnetic 
field^iiSS, or a thin and narrow ring in a perpendicular 
magnetic field^'^. In all cases the transition to the nor- 
mal state is of second order and the vorticity in the wire 
will be equal to zero due to the small cross-section of the 
sample. 

Because the order parameter practically does not de- 
pend on the radial coordinate of the wire we can use the 
one-dimensional model in order to study the response 
of the system on the applied current. In order to take 
into account the suppression of by magnetic field (H) 
we add to the RHS of Eq. (1) the term -{H/Hrytp. 
In some respect it is similar to the way we introduced 
the first type of defect in our wire. We can expect 
that r|^| will increase with increasing H (because the 
'strength' of the 'defect' increases). But because the mag- 
netic field suppresses the order parameter everywhere in 
the sample it also leads to an increase of Ag, because 
in the TDGL model Ag ~ 1/a/|^. It is clear that 
both these processes should decrease jci- The strongest 
mechanism is connected with the change in r|^,|. In- 
deed it is easy to estimate that r|^| ~ 1/(1 — [H/HcY') 
and T0 ~ (1 — (H/Hc)'^)^/'^. Even if we take into ac- 
count that the parameter 7 may decrease with increas- 
ing H (because in non-zero magnetic field there is an- 
other pair-breaking mechanism and instead of te we 
should useiSiia tb/VI + 2tet, with r, = A(T = 0, i? = 
0)/h{H/Hc{T = 0))^ ) it does not lead to an increase 



of jci with an increase of H because T|^| changes faster 
than even in this case. 

In the above model it is easy to show that jc2{H) — 
^4/27(1 - {H/Hc)'^)^-^ for the case of an uniform wire. 
In Fig. 4 we present the results of our numerical calcu- 
lations. We found that both jd and jc2 decreases with 
increasing magnetic field and at some H* they practi- 
cally merge. Unfortunately, it is quite difficult to find an 
analytical expression for the dependence of jci{H) like 
we had for jc2- The reason is the complicated behavior 
of the dynamics of ij: in the phase slip center. 

The variation Ag with increasing H was obtained ex- 
perimentally in Ref.^^. To compare with the theory 
the authors of Ref.^^ used the expressions found in the 
work of Schmid and SchonM, which are valid in the limit 
T ^ Tc- It is interesting that in Ref>i& it was found that 
expressions of Schmid and Schon are quantitatively valid 
even far from T^. From this observation we may hope 
that the present theoretical results are also valid over a 
wider temperature range than only near Tc- 




FIG. 4: Dependence of the lower critical current jd (squares 
and solid curve) and upper critical current jc2 (dots and dot- 
ted curve) on the applied magnetic filed. Dotted curve is 



equation \/4/27(l - 
ted expression 0.26 



(1 - {H/H^rr 



The solid curve is the fit- 
We also plotted the 



dependence jd (H) (open triangles and dashed curve) where 
the effect of // on 7 was taken into consideration. We used 
the simple expression jef / = 7o/ v^l -|- ^o{H/ Hc)^ which even 
overestimates the effect of H (i.e. we took 70 = 'y{H = 0)). 



The main conclusion which follows from the change of 
jd with decreasing wire length and/or applying magnetic 
field is that there exist a critical length L* (or critical field 
H*) below(above) which the current jd becomes equal to 
jc2- It implies that for wires with lengths L < L* and/or 
fields H > H* there will be no jump in the voltage in 
the current-voltage characteristic and the I-V curve will 
be reversible. It will also result in the absence of a S- 
behavior in the V=const regime (see section below). 
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B. Constant voltage regime 

In our earlier work we found that the I-V in the 
V — const regime exhibits an S-Uke behavior— and fur- 
thermore for low voltages there is an oscillatory depen- 
dence of the current on the apphed voltage (see Fig. 5). 
As was shown in Ref— these properties are connected 
with the existence of two critical currents jci and jc2- 
Here, we will discuss the characteristic voltages Vi, V2 
(see Fig. 5) and their dependence on the length of the 
sample. 




0.26 0.28 0.30.,. 0.32 0.34 0.36 




0.31 0.32 0.33 0.34 0.35 




FIG. 5: Theoretical current- voltage characteristics of wires of 
different lengths for 7 = 10. Solid curve corresponds to the 
I=const regime and is practically universal for the considered 
lengths. Dot-dashed curve {L = 40^), dotted curve (L — 80^) 
and dashed curve {L — 160^) correspond to the V = const 
regime, (b) is an enlargement of low voltage region in (a). 

If we apply a voltage V to the wire of length L then 
in the sample an electric field E — V/L exist which will 
accelerate the superconducting electrons. When the cur- 
rent density approaches jc2 phase slip centers will sponta- 
neously appear in the sample. As a result the momentum 
of the superconducting electrons (and hence the current) 
will decrease by 27r/L after each phase slip event. When 
the current density decreases below jd the phase slip pro- 



cess will no longer be active and the applied electric field 
will be able to accelerate the superconducting conden- 
sate. This process leads to periodic oscillations in time 
of the current in the sample. 

We can divide the period of oscillations (at least at 
voltages V < Vi) in two parts. First, the longest part 
is the one during which the condensate is accelerated by 
the electric field till the moment reaches j ~ jc2- The 
second part we call the transition period Ttr- The tran- 
sition period also consists of two parts: the time needed 
for the phase slip processes (which is proportional to the 
number of phase slip events and hence the length of the 
wire) and the time Tq for the decay of the order param- 
eter from •\/2/3 (when j ~ jc2) till the first phase slip 
event and for the recovering of |-0| back to ■\/2/3 from 
zero after the last phase slip event (see Fig. 6). The 
minimal number of phase slip events which occur during 
the transition time is determined by the internal param- 
eters of the superconductor {jd ) and its length^'^ which 
is given by 

N„,i„ = Nint ((pe - Pci)(L/27r + 1)), (9) 

where pd is the smallest real root of the equation jd — 
Pj,i(1— Pel) Nint(x) returns the nearest integer value. 

At y < Vi the time-averaged current {j) increases 
(with oscillations) and in the range Vi < V < V2 \t 
decreases with increasing voltage. The lowest minimal 
current in the latter region is jYi which depends on the 
length of the system (see Fig. 5 and Ref— ). In Fig. 7 
we present the dependence of the above voltages on the 
length of the wire. The explanation for their different 
behavior is the following. At a voltage Vi the period of 
the oscillation T = iirN/Vi (N is the number of phase 
slip events during the transition time) of the current be- 
comes of order 2Ttr- The time Tq does not depend on the 
length and the voltage dX V V\. So, we can estimate 
Vi as 

27riV _ TiN _ TiaL 
' ~ 27^ " To + {tpsc)N ^ To + {Tpsc)aL' ^ 

where (tpsc) is the average time between two phase slip 
events and a is the coefficient which depends on pd (see 
Eq. (9)) and hence on the specific superconductor. When 
To <C {tpsc)<^L the voltage Vi becomes practically inde- 
pendent of the length. Because the time Tq depends on 
the relaxation time of the absolute value of the order pa- 
rameter the length of the wire Lgat at which Vi saturates 
depends on the internal parameters of the superconduc- 
tor. 

The voltage V2 decreases with increasing length of the 
sample because the lower critical current jd decreases. 
It is interesting to note that with accuracy of finding Vi 
the saturated value of the lower critical voltage coincides 
with the voltage jump AV at j — jd in the constant 
current regime. Because the minimal value of V2 is equal 
to Ay for an infinitely long wire we may conclude that 
with increasing wire length the I-V characteristic in the 
range of voltages (Vi, V2) approaches the horizontal line. 
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FIG. 6: The transition time consists of two parts: Tu- ~ 
T0 + T2. 
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FIG. 7: Dependencies of the voltages Vi and V2 on the length 
of the wire. The results are obtained for 7 — 10. 



Defects, magnetic field, short length of the sample etc. 
will decrease the hysteresis and consequently the currents 
jci and jc2 approaches each other. If the difference be- 
tween them will be small enough the S-shape of the I-V 
characteristic at the constant voltage regime changes to 
the usual monotonic behavior (see Fig. 8). This happens 
because when the voltage approaches Vi, the maximal 
current during the period of the current oscillation T will 
be about jd (because jmax cannot substantially exceed 
jc2 but jc2 is already close to jd in this case). 

We also would like to discuss how the change of the 
boundary conditions will change the shape of the I-V 
characteristic. If we apply the N-S boundary conditions 
the TV curve in the voltage driven regime also exhibits 
an S-behavior but without oscillations in the current at 
small voltages (see Fig. 9). In this case there will al- 
ways be an inevitable voltage drop near the boundaries 
connected with the current in the wire through the re- 
lation Vns ^ j^Q (because the electric field and the 
normal current decays on a scale of the charge imbal- 



FIG. 8: Current voltage characteristics of a superconducting 
wire with a local suppression of the critical temperature in the 
center of the wire. The parameters for the defect and the wire 
are the same as in Fig. 2(a). For the TV characteristic with 
small hysteresis in the / = const regime the TV behavior 
in the V = const regime is a single-valued function of the 
current. 



ance length Ag near the boundaries^). Near the bound- 
ary the nonzero electric field is compensated by the term 
AgdE/ds (see Eq. 3(b) - and E = jn in our units) 
and the superconducting electrons are not accelerated 
by this field. The situation is very similar to the case 
when we inject a current in the wire. When the current 
generated by this voltage reaches jci the phase slip pro- 
cess becomes possible in the system. But like the case 
of the current driven regime the superconducting state 
can be stable(metastable) until the current reaches jc2 
at V ~ ic2AQ. From our estimates it follows (see section 
below) that the fluctuations of the order parameter are 
very important in our samples. Therefore, in our calcu- 
lations we introduce at some moment of time (at fixed 
voltage) a phase slip center in the center of the wire and 
checked if it will survive or not. As a result we obtained 
the TV characteristics as presented in Fig. 9. When the 
voltage is less than some critical value {V*) the phase 
slip process decays in time and the resulting current in 
the wire is time- independent. The whole voltage drop 
occurs near the boundaries. At V > V* (and j > jci) 
the dynamics of the condensate in the wire will be similar 
to the case considered above for V > Vi. 

The reason for this is as follows. The current density 
in the wire at small voltages will always be less than jd 
and jc2 due to the relation j ~ Vsn/^q = V/^Q- When 
the current density reaches jd the phase slip process be- 
comes possible in the sample. But such a phase slip pro- 
cess leads to a finite voltage. As a consequence the volt- 
age drop at the boundaries will sharply decrease when a 
phase slip center is created. But it implies that the full 
current density will also decrease and consequently will 
become less than jd- We can conclude that the phase 
slip process may survive in such a type of superconduc- 
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tor only if the applied voltage will be roughly equal to 
AV{j) plus the voltage drop near the boundaries Vsn 
necessary for the creation of a current larger than . 




FIG. 9: Current- voltage characteristics of wires of different 
lengths in the case of the N-S boundary conditions. 

Concluding the theoretical part, we present in Fig. 
10 the I-V characteristics in the I=const and V^const 
regimes at different temperatures close to Tc. With de- 
creasing temperature the range of currents, where there 
is a S-behavior increases and the voltages Vi and V2 in- 
creases. It resembles the experimental results presented 
i n^^'— and in the subsequent Section III. 
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FIG. 10: Current voltage characteristics of a superconducting 
wire of length 600^(0) at different temperatures close to Tc. 
We used typical parameters for Pb: ^(0) ~ 40nm and 7(0) ~ 
100. This result shows an increase of the difference jc2 — jci 
in absolute units with decreasing temperature. 



III. EXPERIMENT 

Measurements of the I-V characteristics were done on 
single Pb nanowires ^^'^'^ with typical diameter about 



50nm and length 22 — (see Table I). At first we 

would like to discuss the rather wide resistive transitions 
and strong magneto-resistance effect (see Figs. 11-12) 
in our samples. We may claim, taking into account the 
small diameter of our nanowires, that the effect of the 
thermo-activited^ and the quantum-activated2SiS& phase 
slip phenomena is very strong in our samples. It is easy 
to see that with increase of the diameter of the nanowire 
the width of the resistive transition decreases. Our esti- 
mations, based on Giordano expressions showed that 
for sample A even at T = the number of quantum phase 
slip events should be about ~ 10^ per second. That is 
the reason why we did not observe in our experiment any 
hysteresis in the current driven regime (see below). But 
this rate is not large enough for destroying the S-behavior 
in the voltage driven regime. Indeed, the period of oscil- 
lations of the current for Pb is about lO^^s a.t V ^ Vi 
and T = 0. It means that only at temperatures close to 
Tc the fluctuating PSC's will interrupt the internal tem- 
poral oscillations in the order parameter and ruin the 
S-behavior. In other words the system can be in a state 
with j > jci only during a time which is less than the 
time between two phase slips. This results in the coinci- 
dence of the I-V characteristics both in the current and 
the voltage driven regimes at temperatures close to Tc 
(see Fig. 3 in Ref,^^). We would like also to mention the 
small difference in the critical temperatures which im- 
plies that all our samples have practically the same value 
of the superconducting gap. 



1.2 




4.5 5.0 5.5 6.0 6.5 7.0 7.5 8.0 8.5 9.0 9.5 

T{K) 



FIG. 11: Resistive transition of our different samples. Note 
that for the narrowest sample (A) the width of the transition 
is widest and for the widest sample (D) it is narrowest. 

As can be seen in Fig. 12, the wider the transition 
width in R{T), the smaller the magnetic fields where 
the resistance starts to increase with increasing magnetic 
field. 

Of course, our nanowires are not free from imperfec- 
tions. Although our fabrication technique produce sam- 
ples which are quite regular (see KeiM^ for Scanning Elec- 
tron Microscopy (SEM) picture), we may distinguish two 
kinds of imperfections. Firstly, there can be deviations 
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TABLE I: Parameters of the different samples. 





L(/im) 


d(nm) 


R„(7.1K)(n) 


R.,,(4.3K)(n) 


He(4.5K)(T) 


5i(4.5K)(nm) 


C2(4.5K)(nm) 


p„{7. IK) {fin ■ cm) 


A 


22 


40 


300.9 


14.9 


1.271 


37 


16 


1.67 


B 


50 


55 


210.2 


21.7 


0.925 


38 


19 


0.86 


C 


50 


55 


465.9 


80.7 


1.652 


21 


14 


1.75 


D 


50 


70 


94.6 


17.6 


0.591 


46 


24 


0.52 




H/H 

c 

FIG. 12: Magneto-resistance of our different samples at 4.3 
K. 

from an ideal cylindrical shape. Indeed, although the 
nanopores of the membrane in which the nanowires are 
grown are designed to be as regular as possible, there may 
exist a smooth variation in diameter along the length of 
the nanowire, which we estimated to be not more than 
10 %. However, in addition to these smooth variations, 
it may happen that a defect in the membrane leads to a 
constriction in the measured nanowire. In this case, the 
diameter can be locally reduced quite importantly and 
this results in a lower critical current and a larger crit- 
ical magnetic field. This is in fact what we observed in 
sample C. In addition to such shape imperfections, we 
note also that the parameters of our nanowires strongly 
vary from sample to sample - see table I (we estimated the 
coherence length using the expression He — 2.9<i>o/(7r^c?) 
(Ci) and He = $o/(27rC^) (6) for the critical field). The 
reason for this difference in resistivity is related to the 
second kind of imperfections, which is structural disor- 
der that is formed during the electrodeposition of these 
nanowires inside the nanopores. Indeed, as shown in 
Ref.^^, these nanowires are polycrystalline and inevitably 
contain structural defects such as dislocations, twins, etc. 
These two kinds of imperfections results in differences in 
the TV characteristics (Fig. 13). For example for sam- 
ples A and B we observed two jumps in the voltage which 
we explain by the appearance of two successive phase slip 
centers in the nanowire but for sample C we found only 
one jump in the voltage. Indeed, due to the presence of 
a constriction in this sample, heating may drastically af- 
fect the behavior of this sample, precipitating the return 



to the normal state. Unfortunately sample D was broken 
during the measurements and we could not measure its 
TV characteristic. 




FIG. 13: Current-voltage characteristics of samples A(T=4.3 
K), B(T=4.37 K) and C(T=4.2K) in the current (grey curves) 
and voltage (black curves) driven regimes. The contact volt- 
age was subtracted from the experimental data. 

We found that jumps in the voltage AV^ are practically 
the same for all three samples (see Figs. 13-14). It agrees 
with the theory presented in Sec. II. Indeed all the three 
samples have almost the same Tc and hence the same 
superconducting gap. It is naturally also to suppose that 
the time te is almost the same for all our samples. We 
can expect that the parameters 7 and u are the same 
for all our samples and consequently the relaxation time 
T|^l and the ratio Aq/^ are the same. This automatically 
leads to the invariance of AV on p„. From the results 
of Sec. lib follows that the voltage V2 decreases with 
increasing (in ^) nanowire length (see Fig. 6). For our 
shortest sample A: T~ 600— 1400^ (see Table I). It means 
that in our nanowires the voltage V2 already reaches the 
minimal value AV^ and hence V2 is independent of p„ for 
our nanowires. The voltage Vi depends not only on the 
length of the nanowire but also on the time change of the 
order parameter during the transition period Tq (see Eq. 
(9)) and hence the voltage V2 may reach Ay at shorter 
lengths than Vi if the time Tq is large enough (in Fig. 6 
the opposite situation is presented - at first Vi reaches 
the saturated value). 

If the current density is uniformly distributed over the 
cross-section of the sample then the oscillations of the 
order parameter will be in phase along the cross-section 
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and in this case the IS.V will not depend on the size of 
the sample. It is a direct consequence of the fact that the 
time-averaged chemical potential of the superconducting 
electrons should be equal and constant on both sides of 
the phase slip center or phase slip line/surface. 

Our results show if one uses the Skocpol-Beasley- 
Tinkham (SBT)^ model for the estimation of Aq one 
should be very careful. Indeed, those authors replaced 
the normal current density in the PS center by the ex- 
pression /„(0) = (/ — (31 ci)- As a result the deriva- 
tive dln(0)/dl — 1 becomes current independent. But 
in general the time derivative may be large than unity 
(see Figs. l(a,b)). Secondly, if the length of the sam- 
ple is comparable with Aq we should replace 2Aq by 
2AQtanh(L/2AQ) (see Eq. (7)). As a result the differ- 
ential resistance in samples with L ~ Aq may be larger, 
in general, than the normal one^S. This case corresponds 
to our samples {Rdif — 320 Ohm, R^if — 351 Ohm and 
Rdif ~ 456 Ohm for samples A,B,C respectively after the 
first jump in voltage). But from the SBT model follows 
that Rdi f < Rn (the equality sign holds for samples with 
L<Aq). 

Unfortunately we do not know the actual dependence 
/„(0)(/) for our samples. If we use the values obtained 
from the SBT model (Aq = 12.6/ito, 48.5/ito and 30.3/zm 
for samples A,B and C, respectively) which is qualita- 
tively understandable for longer samples, we have a too 
small number of PSC (compare samples A and C). And 
in this case, the question arises why Aq changes so much. 
Probably, the SBT model gives us only the correct order 
of magnitude for Aq which is only useful as an estimate. 



only one jump in the voltage exists. We explain it by 
an increase of Aq at relatively large magnetic fields and 
hence there is a lack in space for the coexistence of two 
phase slip centers in the nanowire. At high magnetic 
fields the order parameter is strongly suppressed by H 
and the effect of quantum phase-slip fluctuations be- 
comes more pronounced. This is the reason for a smooth- 
ing of the I-V characteristics at high magnetic fields. 
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FIG. 15: Current-voltage characteristics of sample A (T=4.3 
K) in the current driven regime in the presence of a parallel 
magnetic field. The magnetic field increases from right to left 
with step of O.IT. 
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FIG. 14: Jump in the voltage AV (see Fig. 5(a) and 13) for 
samples A-C at different temperatures. 

Finally, we present our results on the infiuence of an 
applied magnetic field on the I-V characteristics. We 
limit ourselves to data for sample A in the current driven 
regime (see Fig. 15). These results already support our 
theoretical predictions of Section lib. It is evident that 
the lower critical current density decreases with increas- 
ing applied magnetic field. From some range of H- values 



IV. DISCUSSION 

In conventional superconductor near Tc the N-S 
boundary conditions (in the sense mentioned in Sec. 
Ila) are validii. For T — > the bridge boundary 
conditions are more applicable due to Andreev reflec- 
tions at the border between the normal metal and the 
superconductoriiSiSS. At intermediate temperatures we 
expect a mixture of the N-S and bridge boundary con- 
ditions. It means that part of the voltage will drop at 
the boundaries and part in the sample. The situation 
in our experiment is even more complicated because in 
our two-point measurements there is also a voltage drop 
at the contacts. Probably, this will not allow us to ob- 
serve the oscillations of the current in the voltage driven 
regime. Another reason is that our samples are very long 
compared to ^ (even for our 22 fim sample L ~ 1000^) 
and consequently the amplitude of the oscillations is very 
small. But nevertheless our calculations in both limiting 
cases of bridge and N-S boundary conditions predicts an 
S-shape of the I-V characteristic. 

We found that the type of boundary conditions are not 
so important for the process of nucleation of phase slip 
centers if the length of the sample is much larger than 
Aq. However, the difference in the process of the con- 
version of superconducting electrons to normal ones and 
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vice versa at the N-S boundary at various temperatures 
plays a crucial role for the creation of phase slip cen- 
ters in shorter nanowires. It turned out that for similar 
parameters (nanowire's length, coherence length, super- 
conducting gap, te) phase slip centers appear in the su- 
perconducting wire at smaller currents for the case of the 
bridge geometry boundary conditions. 

All our theoretical results are strictly speaking only 
valid near Tc, which is the temperature region where 
Eqs. (1,2) are quantitatively correct. Nevertheless, ex- 
perimental results supports our prediction on the compe- 
tition of two relaxation times in the creation of a phase 
slip center even far from Tc- Indeed, an applied parallel 
magnetic field decreases the lower critical current (com- 
pare Fig. 15 and Fig. 4). Jumps in the voltage AV^ 
turned out to be almost an independent function of the 
disorder (of the resistance of the sample) as it follows 
from theory. Unfortunately, it is rather difficult to check 
the dependence of jd on the length of the nanowire using 
our technique because every preparation of a new sample 
leads to a different level of disorder and hence different 
values for ^ and Aq. 

Therefore, we expect only quantitative differences in 
the dependence ofjci{T, H), jc2{T, H), as compared with 
our theoretical results. Some qualitative differences (see 
Ref^ft) in the dynamic of the order parameter at the 
phase slip center or the creation of charge imbalance 
waves^''' cannot affect the main properties of our theo- 
retical results because it does not influence the existence 
of the two different critical currents: jci and jc2- It may 
lead to quantitative differences in the dependence of 
and T|^| on the microscopic parameters of the supercon- 
ductor. 

Finally, we would like to discuss other mechanisms 
which, for our geometry, may lead to an S-behavior of 
the I-V in the V=const regime. The first is heating2^. In 
order to explain the double S-structure for our samples A, 



B by this mechanism we have to assume that heat dissi- 
pation and heat evacuation have a very complicated and 
non trivial dependence on temperature. We do not know 
any mechanisms which can lead to such adependence in 
our case. Besides we did not any observe hysteresis in the 
current driven regime which is an inevitable property of 
that mechanism if the I-V characteristic would have a S- 
shape in the voltage driven regime. For these reasons we 
believe that heating is not responsible for the observed 
behavior. 

Secondly, if the nanowire contains a normal region the 
I-V characteristic will exhibit an S-behavior due to mul- 
tiple Andreev reflection in the SNS structure'^^. We do 
not have any indication for this process from our R(T) 
and R(H) measurements which shows that our samples 
are quite homogeneous. Furthermore the structure would 
occur at voltages much smaller than A/e. In our case 
the S-behavior is seen for V > A/e - see Fig. 13 (for 
Pb, A(0) ~ lAmeV). Andreev reflection on the N-S 
boundaries may give rise to a zigzag shape of the I-V in 
the V-const regimc'^^ but this effect is negligible for our 
samples with L ^ ^o- 
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